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ABSTRACT 


The increasing number of earth-orbiting satellites, transmitting in 
the 136-138 MHz band, has crowded this allocated frequency region uti- 
lized by the National Aeronautics and Space Administration’s (NASA) 
Space Tracking and Data Acquisition Network (STADAN). STADAN 
ground stations sometimes experience radio frequency interference 
(RFI), caused by the appearance of two or more interfering satellites 
with overlapping 136 MHz emission spectrums, located simultaneously 
within view of a given ground station antenna. 

Satellite orbits can be determined at least one-week ahead; there- 
fore, an a-priori RFI prediction can be made that permits re-scheduling 
of station operations during restrictive periods when RFI is present. A 
computer-controlled satellite RFI prediction technique, employing the 
UNIVAC 1108 digital computer, has been developed that predicts inter- 
ference one-week in advance for up to 50 satellites and 50 stations. 

Computer inputs for the satellites include center frequency, spec- 
tral bandwidth, radiated power level, and Brouwer Mean elements for 
orbit determination. Tracking station inputs include antenna radiation 
patterns, both mainlobe and sidelobes, that are segmented in one-degree 
increments. Satellite-orbit data points, and station- received satellite 
signal-to-interference ratios, are computed at two-minute intervals. 
Station RFI reports show good agreement with predictions. 
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RADIO FREQUENCY INTERFERENCE 
PREDICTION PROGRAM FOR EARTH-ORBITING SATELLITES 

by 

Ralph E. Taylor 

Goddard Space Flight Center 


INTRODUCTION 

The increasing number of earth-orbiting satellites, transmitting in the 136-138 MHz band, has 
crowded this allocated frequency region utilized by the National Aeronautics and Space Adminis- 
tration^ (NASA)' Space Tracking and Data Acquisition Network (STADAN). The 35-odd operational 
satellites, currently supported by STADAN, emit tracking and data acquisition radio signals in the 
136-138 MHz band. The ground stations at times experience radio frequency interference (RFI), 
due to the appearance of two or more interfering satellites having overlapping 136 MHz emission 
spectrums, located simultaneously within view of a given ground station antenna. Department of 
Defense (DOD), other U. S. Government agency and foreign nation satellites, radiating 136 MHz 
signals, also add to the RFI problem. 

A majority of the STADAN station RFI, presently being experienced, is caused by interfer- 
ence from satellites that can be predicted in advance. This is opposed to a minority of random 
electromagnetic compatibility type interference that occurs in an unpredictable manner. 

A control center operator currently employs a satellite visibility chart, that is computer- 
generated, in conjunction with a satellite frequency listing to resolve certain RFI conflicts. This 
report describes an automated computer program for performing essentially the same function 
with increased speed and accuracy. Also, the program calculates the received signal-to- 
interference (S/l) ratio for all visible satellites in addition to making orbit computations. 

Computer inputs include satellite emitted frequency, spectral bandwidth and radiated power 
level. Also inputted are the station geodetic coordinates (longitude, latitude and height), receiving 
system threshold sensitivities, receiving antenna radiation patterns, and Brouwer Mean Satellite 
orbital elements. The program utilizes the UNIVAC 1108 high-speed digital computer written in 
FORTRAN IV format. The UNIVAC 1108 computer was selected since a Brouwer Mean orbit 
generator software program had already been developed for this computer. The high speed and 
relatively large memory capacity (65,000 cores) were also required. 

Certain phases of the UNIVAC 1108 computer program are based upon the deterministic inter- 
ference prediction technique, described by Haber 1 , that calculates the signal-to-interference (S/l) 
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ratio at a given station. Reference 2, on the other hand, determines all visible satellites from orbit 
computation, but does not consider the S/I ratio. 

The RFI prediction program described in this paper can handle only Brouwer Mean orbital 
elements. A future modification of the program will take into account orbits using other than 
Brouwer Mean elements. Also, the program will be expanded for other frequencies including the 
400-401 MHz and 401-402 MHz space bands. 

The RFI prediction program, described in this report, is based upon an initial study effort 
performed by the University of Pennsylvania, The Moore School of Electrical Engineering (Refer- 
ence 1), Philadelphia, Pa.; and a subsequent computer programming effort by the Wolf Research 
and Development Corporation (Reference 3), Applied Sciences Department, Bladensburg, Md. 

DESCRIPTION OF CONCEPT 

Since the satellite orbits can be predicted at least one-week ahead, an RFI prediction can be 
made which makes possible in-advance re-scheduling of station operations during restrictive 
periods when RFI is present. In addition, inferior data taken during an interference interval can 
be identified at the station resulting in data handling cost and manpower savings. 

For a steerable antenna, the satellite being tracked is assumed positioned perfectly within the 
ground station antenna beam (see Figure 1); whereas, the interfering satellite, i, injects an un- 
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STADAN 
STEERABLE 
ANTENNA 
(85-FT.DISH) 

S/I 

TO RECEIVER 


Figure 1 -Radio-link parameters for satellite interference prediction program. 
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desired signal either by transversing the main beam or cutting a high-level sidelobe. The desired 
received signal power, S, and interfering signal power, i, are determined from 
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Range, R 0 , is the independent variable in (1), for a given satellite, assuming an omnidirectional 
transmit antenna gain, G to = 1, constant transmitted power, P to , fixed wavelength, \ 0 , and using 
a steerable ground antenna with a maximum pattern gain, G ro . On the other hand, the received 
interference power level, I , in (2) is affected by the antenna pattern gain variable, G r ., variable 
range, R 4 , and transmitter power, P ti . 

The satellite distance from a given ground station is determined at 2-minute intervals as part 
of the orbit computation for all satellites and ground stations. The appropriate satellite radio-link 
parameters are inserted into (1) and (2) and the station- received signal levels are then calculated. 
These levels are then compared to the threshold sensitivity for the particular receiving system. 
Levels below threshold are automatically eliminated. S/l ratios for all satellites, compared to the 
one being tracked, are then computed from (1) and (2). The magnitude of the S/l ratio, compared 
to a preset value, determines if interference is present or absent. An interference condition is 
defined as existing when the heuristic equality S/l < +20 db holds. 

Prior to computing the S/I ratio, a satellite visibility search sub- routine eliminates satellites 
below a minimum restricted elevation angle (e.g., 10° and below), and a frequency search is made 
based on the criterion, 




By definition, satellites can interfere that meet the (3) criterion. 


( 3 ) 
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The bandwidths, B 0 and B. , in (3) are the emitted spectral bandwidths of the respective desired 
and interfering satellites; whereas, f Q and f . are the respective carrier center frequencies. An 
ideal, uniformly flat, emitted spectral bandwidth is assumed to simplify the computations. The 
emitted spectrum is assumed symmetrically centered about the carrier frequency. The frequency 
search identifies satellite transmitters as potential RFI sources when the (3) equality holds (i.e., 
overlapping spectrums). 

The spectral width of an unmodulated carrier frequency is defined as 30 KHz (i.e., f ± 15 KHz) 
to account for uncertainties in transmitter frequency, Doppler shift, etc. The assumed spectral 
bandwidth for each modulated signal ranges from 60 KHz to 200 KHz depending upon the particular 
modulation employed for a given satellite. These values have been made somewhat larger than the 
actual emitted spectrums to account for the injection of adjacent-channel interference in off -band 
skirt selectivity regions of the ground- station receiver IF amplifier. 

This program utilizes the Brouwer Mean orbital elements to generate orbits at 2-minute inter- 
vals on a magnetic tape. This tape in turn is used in generating RFI predictions for a 24-hour 
period for up to 50 earth-orbiting satellites and 50 world-wide data acquisition stations. The 
stations include Minitrack, 85-ft. dish, 40-ft. dish, VHF range and range rate, SATAN (satellite 
automatic tracking antenna) and 9-element yagi-array antennas. The computer input also includes 
the satellite identification number, radiated power level, transmitter frequency, emitted spectrum 
bandwidth, and ground station antenna radiation pattern data. 

Each STADAN world-wide tracking facility is actually comprised of several stations. The 
Rosman, N. C. facility, for example, has two 85-ft. dishes, a SATAN system and a VHF range 
rate system. RFI is separately predicted for each of the various receiving system types at a given 
site. A considerable saving in computation time, and quantity of RFI print-out data, has resulted 
by developing a standard "Model System" antenna radiation pattern for systems with similar 
antenna be am width and gain characteristics. Such a fT Model System" was developed that represents 
the SATAN, VHF range and range rate, 40-ft. dish and 9-element yagi antennas in STADAN. 


ORBIT ANALYSIS 

Satellite ephemeris calculations 3 in the RFI prediction program are based on the Brouwer 4 
general perturbations method. 

A Brouwer general perturbations orbit generator calculates the time-independent quantities, 
osculating orbital elements, and the inertial rectangular coordinates of a satellite at a given time. 
The zonal harmonics are accounted for, and orbit-drag terms are considered. The computations 
are referenced to the right ascension of Greenwich at January 0.0 of the year of the start of the 
run. 


x s , y s , and z s are the rectangular coordinates of the satellite in inertial space with the x-axis 
pointing toward the vernal equinox, the z-axis pointing toward the North Pole, and the y-axis form- 
ing a right-handed system. After computing the inertial rectangular coordinates of each satellite 


4 



with the Brouwer orbit generator, the earth-fixed rectangular coordinates (x s , Y s and zj of each 
satellite are computed from 

x , = ( x s + y s )“ cos x s • 

Y s = ( x . + y s )* sin • . 


where 


\ s is the satellite’s east longitude , 


arc tan 


- 6 


g ’ 


and 0 g is the right ascension of Greenwich. 

The earth-fixed rectangular coordinates (x ob , Y ob and Z ob ) are computed for each station 
from the geodetic latitude, a>, of the station (antenna); east longitude, y; and height, h, of the station 
above the reference ellipsoid according to, 


X Qb = (h + v) cos $ cos y 

Y ob = (h + v) cos $ sin y 

Zob = [h + V (l - E 2 )] sin « 

Where 


(l - E 2 sin 2 4 >) m 
b 2 


E 2 = 1 - 


2 ' 


and 


a e is the earth's equatorial radius 
b e is the earth's polar radius. 

The local east (E), north (N), and vertical (Z) vectors are computed for each station from, 

E = - sin y i + cos y j + 0 k 

N = - sin $ cos y i - sin $ sin y j + cos 0 k 

A /V A /V 

Z = cos 0 cos y i t cos <$ sin y j + sin $ k . 
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Satellite visibility calculations are made for each time point, for each station, and for each 
satellite. Defining the station- satellite range vector as, 


D i + D j + D k 

X y J z 


where 


i, j and k are the respective x, y and z unit vectors , 

D = X - X . 

x s o b 

D = Y - Y . 

y s o b 

D = Z - Z . . 

z s o b 

The sine of the elevation angle, EL , is computed from, 

R * Z 

,in (EL) = -jj-p 

R • Z is the dot product of the range and vertical vectors and, |r| = (d x 2 + D y 2 + D z 2 ^ * = magnitude 

of range vector. All satellites, below the station's minimum elevation angle of 10°, are ignored. 
Finally, the earth-fixed station- satellite unit vector is 


R 

lil 


For an X - Y axis antenna mount (see Figure 2), the X and Y angles are determined from 



U • E 

= arc tan 

A “ 


u ■ z 

= arc sin 

(U • N) 


arc tan 
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N~ 

0 

■ E 


e - cos (u • z) 

The angle off the North-Zenith plane for the minitrack system, with a stationary antenna 
pattern, is defined as (see Figure 3) 


arc tan 


0 

E " 

_u " 
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Z- ZENITH 



Figure 2-Local coordinate system for steerable 
x-y axis antenna. 


Z 



SATELLITE 


Figure 3— Minitrack antenna fixed local 
coordinate system. 


and the angle off the East-Zenith plane as 


P 


arc tan 


u ■ 

- N 

_u 

• Z_ 


The terms in parentheses are the dot products of the earth-fixed station unit vector, 0 , and 
the respective direction. 

For a steerable antenna with an x-Y mount, the satellite being tracked is assumed positioned 
perfectly within the main beam of the antenna. The angle of the j th interfering satellite, off the 
x antenna axis, is determined from 

j ? 


, J j 

arc tan — — , 

j 

and the angle off the y. axis as 

X 

arc tan -j— . 

j 


COMPUTER PROGRAM 

The UNIVAC 1108 computer program utilizes the following inputs: 

Satellite Parameters (a) Brouwer Mean orbital elements, epoch time, 

and drag terms. 

(b) Transmit carrier frequency (maximum of 3 
per satellite). 

(c) Spectral bandwidth. 
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ANTENNA GAIN - DB 


Ground Station Parameters 


(d) RF power level. 

(a) Geodetic coordinates. 

(b) Antenna radiation patterns in 1 ° increments 
(see Figure 4 example). 

(c) Restricted elevation angles (i.e., 10° and below). 

(d) Receiving system threshold sensitivity. 

(e) S/I criterion. 



E- PLANE PATTERN H- PLANE PATTERN 

Figure 4—136 MHz mini track antenna patterns. 


The above satellite and ground station parameters are entered on punched cards and read into 
the UNIVAC 1108 computer (see Figure 5 flow diagram). Each satellite orbit data point is com- 
puted at 2-minute intervals and the following computer functions then performed in the order listed 
below for each station: 

(1) Satellite visibility search. 

(2) Satellite frequency search. 

(3) Calculation of received signal level and comparison with threshold sensitivity of system. 

(4) Computation of S/l ratio and determination of RFI. 

The UNIVAC 1108 computer then writes the RFI events at 2-minute intervals on magnetic 
tape (see Figure 5), and a station-ordered sort is thereafter performed to provide a summary and 
a detail listing of the predicted RFI events. 

Another magnetic tape is prepared by the UNIVAC 1108 computer, for the Stromberg Carlson 
(SC) type SC 4020 plotter, that graphically plots both a summary and detailed display of the pre- 
dicted RFI events for a single 24-hour period. The same data is available on both 35 mm film (for 
optical viewer) and 9-inch wide rolled-paper chart printouts. 

The UNIVAC 1108 also writes a listing of the predicted RFI events on magnetic tape for the 
IBM 1401 or 360 type computer-printer. The latter equipment prints-out a tabulated (15 inches 
wide) IBM listing of the predicted RFI events identifying the satellites involved and their respective 
frequencies. 
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STATION DATA: 


SATELLITE DATA: 



Figure 5— Flow diagram, satellite interference prediction computer program. 


All three data forms (i.e., 35 mm film, SC 4020 chart-roll and IBM paper printout listing) are 
available as a data set, covering a given 24-hour period, for use in controlling STADAN scheduling. 
Seven (7) RFI prediction data sets (24-hrs. per set) are required to cover a seven (7) day week. A 
week's worth of RFI data can be run at one time. The total UNIVAC 1108 computer running time 
for the RFI prediction program is about one hour per week. This amount of computer time is 
not considered excessive for operational use. The RFI prediction program output data (SC 4020 
plots and IBM 1401 printout) can be used by Operations Control Center personnel for determining 
clear-channel station scheduling. 
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DATA PRINT-OUT FORMAT 


SC 4020 Plotter 

The SC 4020 plotter output includes a summary bar chart (see Figure 6 example) which lists 
sequentially-numbered ground stations on the ordinate, and Zulu time on the abscissa in 1-hour 
increments from 0 to 24 hours. Each predicted RFI event is indicated by the printed symbol f T’, 
and a solid black-bar section (made by repeatedly printing symbol T T") indicates continuous inter- 
ference for that time period. The summary bar chart is useful for an initial tT quick look" check to 
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Figure 6— SC 4020 plotter summary of predicted satellite interference. 
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identify RFI regions in STADAN. Figure 6 gives the WORST CASE condition since the program 
assumes that all satellites continuously radiate. The indicated RF1, for a given station, represent 
contributions from all satellites in the program. 

Detailed graphs are also available on the SC 4020 plotter printout that lists the abbreviated 
names of the satellites on the ordinate (sequentially numbered) and time from 0 to 24 hours on the 
abscissa (see Figure 7). The satellites involved in a particular RFI event are identified by a 
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Figure 7— SC 4020 plotter detail display of predicted satellite interference. 
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26 ROSMAN 85 FOOT DISH MAY 22, 1967 


vertically printed arabic number,* located opposite the satellite name abbreviation, at a particular 
time. A solid-bar section again indicates continuous RFI for that region. 

The horizontal thin-line marked regions, sometimes appearing as dots on the detailed SC 4020 
graph (Figure 7), indicate the period of time a particular satellite is visible over a given station. 

A satellite, by definition, must always be above the horizon to contribute in an RFI event at a given 
station. However, two or more satellites, simultaneously visible over a given station, may or may 
not contribute to an RFI event depending upon their frequency, and the resulting signal-to- 
interference (S/l) ratio. Using a station- scheduling criterion that eliminates passes for all visible 
satellites, that are close in frequency without regard for the S/l ratio, unnecessarily eliminates a 
significant amount of good satellite signal data. 

The detailed SC 4020 graphs are sequentially listed wherein generic equipment types are 
grouped together for convenience in viewing (e.g., all Minitrack systems in sequence, 85-foot 
dish, etc.). 


IBM 1401 or 360 Printout 

The IBM 1401 or 360 Tabulated printout lists satellite input quantities to the computer (see 
Figure 8), input parameters for the tracking stations (Figure 9), and a typical printed-page listing 
of the predicted RFI events (see Figure 10). 

SATELLITE INPUT QUANTITIES 


BROUWER MEAN ORBITAL ELEMENTS TRANSMITTING 

CHARACTERISTICS 


SAT 
IDENT . 

EPOCH 

YYMMOO HHMMSS 

SEMI-MAJOR ECCEN 
AXIS (Ml 

INCL 

MEAN 

ANOMALY 

ARG OF 
PERIGEE 

RA OF 
NODE 

POWER 

WATTS 

FREQ BNO WTH 
(MHZ) (KHZ) 

PEG-A 

670426 

OCCOOO 

6988339 .015782 

31.766 

315.032 

133.153 

289.301 

M 

136.889 

136.410 

60. 

120 . 

PEG-B 

670423 

cccooc 

6997593 .015376 

31.776 

219.647 

229.682 

198.114 

• 100 
.750 

136. 889 
136.410 

60 . 

120 . 

PGG-C 

670419 

CCQQ0C 

6892 794 .001622 

28.889 

61.607 

180.235 

11.633 

• 100 
.750 

136.592 
136.4 10 

60. 

120 . 

RELAY2 

670424 

CCCOOC 

11129457 • 24C8G7 

46.349 

155.296 

59.325 

350.879 

• 250 

• 250 

136.621 

136.621 

30. 

120 . 

SANMAR 

670426 

101154 

6862142 .038605 

2.892 

348.184 

296.291 

131.785 

• 100 

136.740 

30. 

SYNC-2 

670315 

153439 

42135819 .000386 

30*356 

359.312 

.639 

298.429 

2.000 

2.000 

2.000 

136.470 

136.47C 

136.980 

30. 

120 . 

120 . 

SYNC- 3 

670315 

155733 

42161964 .000185 

1.530 

253.598 

286.535 

91.253 

2.000 

2.000 

2.000 

136.470 

man 

30. 

i!8: 

TIR0S7 

670415 

CCCOOO 

7011365 .002122 

58.232 

193.019 

123.529 

107.082 

• 050 

• 050 

136.233 

136.921 

120 . 

120 . 

TIRQS 8 

670408 

CCOOOC 

7104856 .003455 

58.497 

86.115 

173.884 

162.191 

• 050 
.050 
5.000 

136.23 2 
136.924 
136.924 

120 . 

118 : 

TIR0S9 

670425 

ocooco 

8018437 .116932 

96.418 

10.432 

214.054 

39.100 

:8io 

136.232 

136.919 

118 : 
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670425 

CCOOOC 

7164599 .006731 

98.579 

126.283 

106.614 

16.605 

• 050 

• 050 

136.232 

136.925 

120 . 

120 . 

UK- 2 

670513 

142912 

6902538 .037001 

51.639 

183.950 

176.413 

23.093 

.250 

136.560 

90. 

UK- 3 
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16C932 

6925857 .007641 

80.182 

352.600 

172.108 

170.082 

.250 

136.560 

90. 


Figure 8— Satellite input quantities to UNIVAC 1108 computer. 

_ 

*e.g. 3 for 13 . 
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TRACKING STATIONS 



NAME 

TYPE 

GEODETIC LAT 

EAST LONG 

HE IGHT 

MIN 

THRESH 

GAIN 

S/I 





DO 

MM 

SS 

DD 

MM 

SS 

I METERS) 

EL. 

(DBM) 

(DB) 

<D8) 

1 

FTMYRS 

MINITRACK 

26 

32 

5 3*52 

278 

8 

3.89 

3.7 

10.0 

-132.0 

16.30 

60.0 

2 

JOBURG 

MINITRACK 

-25 

52 

58.86 

27 

42 

27.93 

1521.0 

10.0 

-132.0 

16.30 

60.0 

3 

LIMAPU 

MINITRACK 

-11 

46 

36.49 

2 82 

50 

58.18 

49.1 

10.0 

-132.0 

16.30 

60.0 

4 

MADGAR 

MINI TRACK 
M IN ITRACK 

-19 

0 

27.11 
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Figure 9-STADAN tracking station input quantities to UN1VAC 1108 computer. 
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Figure 10-Typical IBM 1401 printout listing of predicted satellite interference. 
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The SC 4020 plotter output can be initially checked to determine "quick-look" RFI conditions 
in STADAN for a given projected 24-hr. period; whereas, the IBM printed-page (see Figure 10 
example) can be consulted for an accurate time listing of the predicted RFI events. This listing 
includes Zulu Time (hours, minutes, seconds - HHMMSS), satellite abbreviated name, and corre- 
sponding carrier center frequency(s). 

The existing UNIVAC 1108 computer RFI prediction program can handle up to a maximum 
of eight (8) RFI events at any one instant of time; therefore, there is an upper limit of eight 
(8) RFI events that can be printed for any given minute of time. Occasionally, the IBM print- 
out lists two seemingly redundant RFI events during the same minute. This indicates two sep- 
arate RFI events caused by broadband modulation, or a change in signal power level, for the 
same center frequency. 

SATELLITE RFI PREDICTION TEST RESULTS 

In order to check the program prediction accuracy, several trial-run RFI predictions have 
been made including both STADAN- supported and non- supported satellites. The initial computer 
run, made on March 29, 1967 with 33 satellites, verified four (4) out of eleven (11) station- 
reported RFI events involving satellites. In one instance, Santiago commanded-on ALOU-1, 
radiating on 136.077 MHz, while simultaneously tracking IMP-3 on 136.124 MHz. The two signals 
were close enough in frequency to cause interference, and the predictions revealed that ALOU-1 
would interfere with IMP-3 from 29/0355Z to 0405Z. Santiago reported that RFI existed from 
29/0352Z to 0400Z and further suspected ALOU-1 as the interfering source. This prediction 
program would eliminate future RFI events of this type by re-scheduling of passes. 

A later prediction run, made on May 25, 1967 with 43 satellites, had an improved score 
wherein eight (8) out of eleven (11) satellite RFI reports were verified. A still later opera- 
tional run, with participating Operations Center Branch (code 512) personnel, was made cover- 
ing the period September 19-25, 1967 that included 40 satellites and 27 STADAN network sta- 
tions. This run was more successful and resulted in a better score wherein 40 out of a 
possible 47, or 85%, of the station- reported RFI events were accurately predicted (see Table 
1). The interference time intervals, as reported by the stations, agrees closely with the pre- 
dicted time intervals. In several instances, the predictions reveal that three (3) satellites 
could have been involved in a single RFI event. 

A rather large amount of potential interference is indicated by the SC 4020 summary 
chart (see Figure 6 typical excerpt). The RFI picture improves by observing a typical de- 
tailed display output, from the SC 4020 plotter (see Figure 7), that shows predicted interfer- 
ence for the 85-foot dish antenna at Rosman, N. C. The predicted interference is somewhat 
larger than actually exists in STADAN since the computer program assumes that all 136 MHz 
emitters are continuously radiating whether commandable or not. 
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Table 1 

Forty (40) Predictions Verified by ST AD AN RFI Reports for September 19-25, 1967. 


RFI Prediction 

Station- Reported Interference Event 

IQSY-1 (136.530) with ERS-18 (136.530) at 19/0118Z to 0120Z 
ATS- 2 (136.470) with ERS-18 (136.530) at 19/0758Z to 0940Z 
RELAY-1 (136.140) with IMP-4 (136.140 at 19/0814Z to 0824Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 19/0444Z to 0550Z 
NIMB-2 (136.497) with ERS-18 (136.530) at 19/1150Z to 1158Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 19/1100Z to 1144Z 
ATS-2 (136.470) with ERS-18 (136.530) at 19/1158Z to 1236Z 
ESSA-5 (136.770) with EGRS-7 (136.800) at 19/1710Z to 1720Z 

FR-1A (136.800) with ESSA-2 (136.769) at 20/0054Z to 0102Z 
ATS-2 (136.470) with ERS-18 (136.530) at 20/1130Z to 1304Z 
DME-1 (136.380) with ERS-27 (136.380) at 20/1342Z to 1350Z 
NIMB-2 (136.498) with ERS-18 (136.530) at 20/0400Z to 0408Z 
ESSA-4 (136.770) with ESSA-2 (136.769) at 20/0506Z to 0516Z 
PEG-C (136.591) with RELAY-2 (136.621) at 20/0846Z to 0850Z 

UK-3 (136.560) with IQSY-1 (136.530) at 21/0916Z to 0918Z 
IQSY-1 (136.530) with UK-3 (136.560) at 21/1056Z 
ERS-18 (136.530) with IQSY-1 (136.530) at 21/1056Z \ 

UK-3 (136.560) with IQSY-1 (136.530) at 21/1056Z J 

RELAY-1 (136.140) with IMP-4 (136.140) at 21/0622Z to 0714Z 

DME-1 (136.380) with ERS-27 (136.380) at 22/0412Z to 0416Z 
ALOU-1 (136.592) with ERS-18 (136.530) at 22/1606Z to 1608Z 
IQSY-1 (136.530) with ERS-18 (136.530) at 22/2316Z to 2322Z 1 

ATS-2 (136.470) with ERS-18 (136.530) at 22/2322Z to 2326Z J 

RELAY-1 (136.140) with IMP-4 (136.140) at 23/1412Z to 23/1414Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 23/1734Z to 1754Z 
UK-3 (136.560) with ERS-18 (136.530) at 23/1948Z to 1954Z V 

ALOU-2 (136.587) with ERS-18 (136.530) at 23/1956Z to 2004Z J 
ESSA-1 (136.918) with TIROS-7 (136.921, at 23/1702Z 

FR-1A (136.800) with ESSA-5 (136.770) at 24/0108Z to 0112Z 
IQSY-1 (136.530) with ERS-18 (136.530) at 24/0740Z to 0744Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 24/1426Z to 1520Z 
IQSY-1 (136.530) with ERS-18 (136.530) at 24/0922Z to 0930Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 24/0506Z to 0608Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 24/0830Z to 0842Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 24/2040Z to 2116Z 

FR-1A (136.800) with ESSA-2 (136.769) at 25/0310Z to 0316Z 
ALOU-1 (136.592) with UK-3 (136.560) at 25/0852Z to 0856Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 25/1134Z to 1248Z 
IQSY-1 (136.530) with ERS-18 (136.530) at 25/23 56 Z 
IQSY-1 (136.530) with ERS-18 (136.530) at 25/2212Z to 2218Z 
IQSY-1 (136.530) with ERS-18 (136.530) at 25/1918Z to 1926Z 
RELAY-1 (136.140) with IMP-4 (136.140) at 25/0554Z to 0650Z 
ATS-2 (136.470) with ERS-18 (136.530) at 25/0850Z to 0914Z 

IQSY-1 with ERS-18 at 19/0119Z (JOBURG) 

Satellite RFI with ERS-18 at 19/0830Z to 0834Z (JOBURG) 
RELAY- 1 with IMP-4 at 19/0816Z (JOBURG) 

RELAY-1 with IMP-4 at 19/0500Z to 0521Z (MADGAR) 

Satellite RFI with ERS-18 at 19/1151Z to 1159Z (JOBURG) 
RELAY- 1 with IMP-4 at 19/1 100Z to 1130Z (SANTIAGO) 

Satellite RFI with ERS-18 at 19/1226Z to 1230Z (JOBURG) 
Satellite RFI with EGRS-7 at 19/1 709Z to end of pass (ALASKA) 

Satellite RFI with ESSA-2 at 20/0100Z (LIMA PU) 

Satellite RFI with ERS-18 at 20/1236Z to 1238Z (JOBURG) 

DME-1 with ERS-27 at 20/1342Z to 1347Z (JOBURG) 

NIMB-2 with ERS-18 at 20/0359Z to 0405Z (ORRORAL) 

Satellite RFI with ESSA-2 at 20/0508Z (ALASKA) 

RFI with RELAY-2 at 20/0844Z to 0845Z (SANTIAGO) 

UK-3 with IQSY-1 at 21/0917Z (JOBURG) 

IQSY-1 with UK-3 at 21/1054Z to 1056Z (JOBURG) 

UK-3 with IQSY-1 at 21/1050Z to 1059Z (JOBURG) 

RFI with IMP-4 at 2 1/061 5Z to 0628Z (MADGAR) 

Satellite RFI with ERS-27 at 22/0412Z to 0417Z (ORRORAL) 
Satellite RFI with ERS-18 at 22/1605Z to 1609Z (FT MYRS) 

IQSY-1 with ERS-18 at 22/2318Z to 2323Z (JOBURG) 

RFI with IMP-4 at 23/1407Z to 1414Z (FT MYRS) 

RFI with IMP-4 at 23/1736Z to 1748Z (FT MYRS) 

Satellite RFI (UK-3) with ERS-18 at 23/1951Z to 1956Z (ALASKA) 
RFI with TIROS-7 at 23/1701Z to 1703Z (SANTIAGO) 

Satellite RFI with ESSA-5 at 24/0105Z to 0114Z (ALASKA) 

IQSY-1 with ERS-18 at 24/0742Z to 0746Z (JOBURG) 

RELAY- 1 with IMP-4 at 24/1451Z to 1510Z (SANTIAGO) 

IQSY-1 with ERS-18 at 24/0922Z to 0932Z (JOBURG) 

RELAY- 1 with IMP-4 at 24/0520Z to 0530Z (MADGAR) 

RELAY- 1 with IMP-4 at 24/0830Z to 0837Z (MADGAR) 

RELAY- 1 with IMP-4 at 24/2040Z to 2045Z (ORRORAL) 

RFI with ESSA-2 at 25/03 10Z to 0314Z (MADGAR) 

Satellite RFI with UK-3 at 25/0855Z (JOBURG) 

RELAY-1 with IMP-4 at 25/1153Z to 1239Z (SANTIAGO) 

IQSY-1 with ERS-18 at 25/23 54Z to 26/0004Z (ORRORAL) 

IQSY-1 with ERS-18 at 25/2211Z to 2219Z (ORRORAL) 

IQSY-1 with ERS-18 at 25/1917Z to 25/1923Z (ROSMAN) 

RELAY- 1 with IMP-4 at 25/0610Z to 0641Z (MADGAR) 

Satellite RFI with ERS-18 at 25/0852Z to 0857Z (SANTIAGO) 
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The space above the earth is becoming saturated with 136 MHz emissions as revealed from a 
prediction of the number of multiple satellites, interfering simultaneously, over a typical STADAN 
station. Twenty-nine (29) near-earth satellites were considered using the MADGAR (Tananarive, 
Malagasy Republic) Minitrack site. The RFI prediction printout listing for MADGAR revealed a 
total of 103 potential RFI events, involving two (2) or more satellites, for the twenty-four (24) hour 
period on March 29, 1967 

There was one unusual RFI event predicted for the MADGAR Minitrack station involving four 
(4) satellites, simultaneously (see Table 2). If the MADGAR site had been tracking the BE-B satel- 
lite on 136.171 MHz, four (4) minutes of continuous interference from three (3) other satellites 
could have existed from 29/OOOOZ to 29/0004Z. 


Table 2 

Predicted RFI Events Involving Four (4) Satellites for Madgar Minitrack 
Tracking BE-B Satellite. 


RFI 

Event Time 
(HH MMSS) 

Station: Madgar Type: Minitrack Date: March 29, 1967 

Satellites Involved in Single RFI Event 

29/000000 Z to 000400 Z 
(4 minutes) 

OGO-l (136.201 MHz ± 15 KHz) with BE-B (136.171 MHz ± 30 KHz) 
OGO-3 (136.200 MHz ± 15 KHz) with BE-B (136.171 MHz ± 30 KHz) 
IMP-3 (136.124 MHz ± 20 KHz) with BE-B (136.171 MHz ± 30 KHz) 


STADAN-STATION YEARLY RFI REPORT SUMMARIES 

Yearly summaries of STADAN- station RFI reports, for the past three (3) years, indicate a 
high percentage of the reported interference was due to RFI from satellites (see Figures 11, 12 
and 13). For example, the Table 3 summary shows that from 58% to 94% of the reported RFI 
events were caused by satellites for the three (3) stations considered. 

The average time duration of the station-reported RFI events in 1966 ranges from 0.41 minutes 
per event for the Lima, Peru, station to 11.9 minutes per event for the Fort Myers, Florida, station. 

A word should be said at this point about the total number of reported satellite RFI events, in 
the yearly summaries, that range from seven (7) to 115 interference events per station, in 1966, 
for the three (3) stations considered. The computer-predicted interference for the MADGAR 
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NO. INCIDENCE PER YEAR MO. INCIDENCE PER YEAR 




YEAR 

AVERAGE DURATION OF INTERFERENCE 
( minutes) 

YEAR 

AVERAGE DURATION OF INTERFERENCE 
(minutes) 

1964 


1964 


1965 

23.1 

1965 


1966 

11.9 

1966 

0.41 


TOTAL AVERAGE 12.27 

Figure 11— Yearly RFI report summary (Ft. Myers, 
Florida station). 



AVERAGE DURATION OF INTERFERENCE 
( minutes) 


1964 6.25 

1965 3.4 

1966 7.08 

TOTAL AVERAGE 6.59 


Figure 12— Yearly RFI report summary (College, 
Alaska station). 


TOTAL AVERAGE 0.41 

Figure 13— Yearly RFI report summary (Lima, 
Peru station). 


Table 3 

Summary of STADAN-Station RFI Events for 1966. 


STADAN 

Station 

Total No. 
RFI Events 

No. Satellite 
RFI Events 

Percentage 
RFI Events 
Due to Satellites 

Fort 
Myers , 
Florida 

29 

21 

72% 

College, 

Alaska 

122 

115 

94% 

Lima, 

Peru 

12 

7 

58% 


Minitrack station, for example, indicates a much 
larger number of RFI events (i.e., total of 103 
events per twenty-four (24) hours). There are 
several reasons for the apparent difference; 
namely: 

(1) For analytical reasons, the computer- 
prediction program assumes that all 136 
MHz satellite transmitters are contin- 
uously radiating. 
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(2) Satellite "on-off” command schedules reduce the satellite transmitter "on" time. 

(3) STADAN stations do not track or take data on all satellite passes. 

The number of STADAN- station radio frequency interference reports have substantially in- 
creased in 1967 compared to 1966. For example, during the month of May 1967, 358 separate 
interference reports were filed by the network stations. The largest portion, 271, was identified 
or suspected as being caused by multiple satellites, radiating in the 136-138 MHz band, simultane- 
ously visible over the same site. This figure is much larger than the value of only 68 satellite 
interference reports, from all stations, for the month of December, 1966. The remaining 87 inter- 
ference events for May 1967 were identified as being caused by solar flare activity, aircraft, or 
from unknown sources. 


CONCLUSIONS 

A computer-controlled prediction program has been developed for the STADAN network to 
reduce radio frequency interference between multiple satellites by means of appropriate station 
operations re- scheduling. This program, written in Fortran IV format for the UNIVAC 1108 
operational computer, will handle up to 50 satellites and 50 ground stations. Forty (40) satellites, 
including certain non-supported satellites, have been included in the program to date. 

A new RFI prediction program is currently being written, in FORTRAN G and H format for 
the IBM 360/75 computer, that will handle up to 100 satellites; accept Brouwer Mean and NORAD 
type orbital elements with equal accuracy; and have as an input the satellite command schedule. 

The existing program is flexible in that predictions can be made for certain selected 
stations and satellites. For example, predictions for the 10 Minitrack stations may be run as a 
separate program; and interference conflict predictions, for as few as 2 satellites, can be made 
for a given 24 hour period. Furthermore, program parameters, including the threshold signal-to- 
interference ratio and satellite transmit spectral bandwidth, can be readily adjusted to meet various 
mission requirements. 

The accuracy of the program has been verified by comparing RFI predictions with STADAN- 
station reported satellite interference, and good agreement has been obtained. In addition, this 
program can be utilized to generate satellite orbit predictions for the STADAN stations. For ex- 
ample, generated orbit prediction printouts were compared with conventional operations-type pre- 
dictions for the GEOS- A, IQSY, ESSA-3, and Relay-2 satellites. Good agreement was obtained 
wherein the antenna x-y axis pointing angles correlated within one-half degree (0.5°) for the two 
methods. 

The highly-elliptical orbit satellites, such as AIMP-D, are not included in the existing com- 
puter prediction program. However, a number of the satellites reported through the North 
American Air Defense Command (NORAD) have already been integrated into the program. These 
include NORAD satellites, UK-2, EGRS-6, EGRS-8, SYNCOM-2 and SYNCOM-3. 
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This prediction program should be effective in reducing the steadily increasing number of 
STADAN RFI events involving multiple satellites. For example, in May 1967, there were about 
270 reported satellite RFI events compared to only 68 satellite RFI events in December 1966. 

The May 1967 data indicates that serious RFI from satellites exists in about 2% of all passes, and 
is steadily increasing as more satellites go into orbit. 
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